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ABSTRACT 



We study the thermal evolution of a pulsar after a glitch in which the energy 
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{T) • is released from a relative compact region. A set of relativistic thermal transport 
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making the assumption of spherical symmetry. We use an exact cooling model to 
solve this set of differential equtions. Our results differ significantly from those 
obtained under the assumption of spherical symmetry. Even for young pulsars 
with a hot core like the Vela pulsar, we find that a detectable hot spot can be 
observed after a glitch. The results suggest that the intensity variation and the 
relative phases of hard X-ray emissions in different epoches can provide important 
information on the equation of state. 

Subject headings: dense matter - stars: evolution - star: interiors - stars: neutron 
- stars: X-rays 



1. Introduction 



The studies of the thermal evolution of pulsars are believed to provide vital informa- 
tion on the internal properties of the neutron stars (for a review see e.g., Tsuruta 1992). 
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Theoretical cooling curves are often compared with the surface temperatures of pulsars with 
different ages. However, different pulsars clearly have different parameters, e.g. the rotation 
period (P) and the surface magnetic field (B), which may affect the cooling processes. Sev- 
eral proposed internal heating mechanisms, e.g. the fractional heating between the crustal 
superfluid and the crust (Alpar et al 1984a, Shibazaki and Lamb 1988), the crust cracking 
(Cheng et al 1992), the chemical heating (Reisennegger 1995) etc., all depend on the pulsar 
parameters, i.e. P and B. Furthermore, different pulsars have different internal properties, 
e.g. mass, equation of state, impurity content etc., which strongly affect the cooling curve. 
Hence ideally speaking, one would like to study the thermal evolution process of an indi- 
vidual pulsar. But realistically, the normal cooling process is an extremly slow process and 
cannot be followed observationally for an individual star. 

The thermal evolution of a pulsar after a glitch is particularly interesting in that it fills 
in this gap. A large amount of energy could be suddenly released during the period jump 
(glitch). The energy will eventually be transported to the surface of the star, and released 
in the form of a transient thermal X-ray emission. Several authors (Van Riper et al 1991, 
Chong and Cheng 1994, Hirano et al 1997) have calculated the thermal evolutions of pulsars 
after glitches, but they all assumed that the energy is generated in a spherical shell inside 
the star. They found that the glitches cannot produce very significant observed results for 
young pulsars like the Vela pulsar (Van Riper et al 1991). However, the energy released in 
the glitch is likely to be deposited in a compact region in the inner crust, by either superfluid 
unpinning (e.g. Cheng et al 1988, Alpar and Pines 1995) or crust cracking (e.g. Ruderman 
1991, Cheng et al 1992). The energy transport is clearly not spherical symmetric (although 
the background geometry is to a good approximation spherical). In this paper, we study the 
general case of the energy transport inside neutron stars after a glitch without the spherical 
symmetry assumption, and arrive at a conclusion substantially different from previous re- 
sults. In section 2, we derive the general relativistic non-spherical symmetric transport and 
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energy balance equations. The necessary physics inputs and solution algarithms for solving 
this set of relativistic differential equations are decribed in section 3. Numerical results are 
presented in section 4, with a brief discussion in section 5. 



2. General Relativistic Non-Spherical Symmetric Thermal Transport and 

Energy Balance Equations 

The Newtonian thermal transport and energy balance equations are given respectively 

by 

V • F = -C v ^ = -nT^ - Q„ (1) 



and 



VT = -1 (2) 



where F is the energy flux, Cy is the specific heat capacity, T is the temperature, n is the 
particle number density, s is the specific entropy per particle, Q v is the neutrino emissivity 
per unit volume and K is the thermal conductivity. 

To generalize the above equations to relativistic situation, we make the following as- 
sumptions. (1) The neutron star is rotating slowly enough that the metric tensor describing 
the background spacetime can be written as (Tolman 1934) 

fa* = _ e 2Z(r) dt 2 + e 2K(r) dr 2 + ^2 + ^2^2) ? ( g) 

where e A ^ = (1 — 2Gm/rc 2 )~ 1 / 2 . (2) The diffusion limit is appropriate for the energy 
transport. (3) There are no other entropy-generating mechanism besides diffusion, and 
second order flux terms in the transport are neligible. (4) Fluid motion inside the star due 
to the thermal effects is negligible. These simplifying assumptions are clearly justified for 
the problem at hand. The energy transport is then governed by the following equations. 
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The energy-momentum tensor can be written as 

T" v = {p + P)UW + pg^ + qW + U^q" (4) 

where If 1 is the 4-velocity of fluid flow, p is total energy density measured in the rest 
frame of fluid, and P is the presure in rest frame of fluid. The heat flow is given by g M = 
-K{g» v + U^U U ){T, V +Ta u ) where a v = U^U^ is the 4-acceleration, and K = K(T) is the 
effective conductivity, e.g. K oc T 3 (t, r, 9, <p) for photon diffusion. Since we assume that there 
is no fluid motion the 4-velocity of the fluid is *7 a =(e _ *,0,0,0) and a a =(0,e" 2A 9 r <l>,0,0). The 
4-heat flow is given by q a = (0, —Ke~ 2A (d r T + T«9 r $), -$d e T, -^§^d v T). The energy- 
momentum tensor components related to the energy transport equation are following: T tt = 
pU'U^pe' 2 ^, T tr = f r = q r U t =—Ke~ ib ~ 2K {d r T + T<9 r $), T te = f = ^U*=-E^.d e T 
and T 1 ^ — f v — q Lp U t = — ^in 2 e ®<pT. To express the thermal transport equation in terms of 
local observables, we note that T a @ = f^f^g-^ ^ where d( a = (e' s 'dt,e A dr,rd8,rsin8d(p) are 
the unit one forms, and dx a =(dt,dr,d8,dip). The locally measured energy flux in spherical 
coordinate is 

F = (T* f , T {§ , T&) (5) 



where 



and 



,t> = e (A+<K) / r = -Ke-*- A d r (e*T), (6) 

T ie = re * f e = _^ df)T ( 7 ) 
r 



= e^rsinOf* = — d v T. 

rsmO 



Therefore, the relativistic thermal transport equation is given by 

(e*T), = (9) 

where ; i denotes spatial covariant derivative on the constant time slice of the metric given 
by Eq.(|3p, with i — r,9, (p. We can see that this equation (i) in the Newtonian limit reduces 
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to to Eq.(0) with $ and A go to zero, and (ii) in the spherical symmetric case reduces to the 
equations of, e.g., Straumann (1984), with J| = j- = 0. 

The energy balance equation can be derived by the conservation of the energy-momentum 
tensor, namely = T\ = T% + + f + d r f r + f r d r A + 3fd r $ + d e f 9 + d v f* 
where T tt . t = ffe -2 * = §fe~' 1 '. Here, |^ is the rate of change of energy density measured 
in proper frame. % depends on the processes under consideration, e.g., if only heat con- 
duction is considered, one has = CVff, with Cy being the heat capacity measured in 
the proper frame. We include neutrino emission, which leads to |j = + e®Q v = 

^^d r {r 2 f r e k+ ^) + j^d e (sin6f d ) + e^d^f*. Using the relations between f a and F a in 
Eqs (§) to (§), we obtain 

dT $ e -(*+ A ) Q 2 . 2$ e* d § e* d 

^ V dt ^ 6 r 2 dr^ e ) + rs i n Q qq^ s ' rs { n Q ^ ^ 

It is very easy to see that the above equation again recovers the well-known spherical sym- 
metric well as the Newtonian limit. 



3. Physics Inputs and Numerical Algarithms 

Physics inputs include stellar model, thermal conductivity, heat capacity, neutrino emis- 
sivity, superfluidity, surface temperature treatment, and the position and amount of the en- 
ergy released by a glitch. In our calculations, we divide a neutron star into two parts. The 
region from the neutron star center to where the mass density equals the nuclear density 
Pn = 2.8 x I0 u g/cm 3 is defined as core, while the region from p^ to the boundary density 
Pb = 10 9 g/cm 3 is defined as crust. The core is treated as isothermal with a time- dependent 
temperature. In the crust region, the temperature after a glitch has both spatial and tempo- 
ral variations; hence an 'exact' or 'evolutionary' treatment is necessary (Nomoto and Tsuruta 
1987, Van Riper 1991). In our treatment, the energy flow and the spatial dependence of the 
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temperature throughout the crust is followed but the thermal effects on the stellar structures 
are neglected, i.e. the same hydrostatic stellar model is used throughout the evolution. 

The stellar structure of a neutron star is determinded by the equation of state (EOS). 
We consider three representative EOSs in this paper, with the total mass of the neutron star 
in all three models taken to be 1.4M (M is the solar mass). The first EOS is BPS (Baym, 
Pethick, and Sutherland 1971), which is often used as a soft EOS in neutron star cooling 
studies. BPS model results in high central density and little mass in crust, and hence a 
small radius. The second EOS used is PPS (Pandharipande, Pines, and Smith 1976), which 
is a stiff EOS with low central density, thick crust and large radius. The third EOS UT 
(Wiringa and Fiks 1988), is a representative intermediate stiff model. Its central density, 
crust thickness and radius are between those of BPS and PPS model. 

In the crust region, we fitted the thermal conductivity data provided by Van Riper 
(1991), who followed the work of Itoh et al. (1984a,b,c,d), together with the quantum 
corrections of Mitake, Ichimaru, and Itoh (1984). The ions A and Z which required in the 
conductivity formulae are taken from Lattimer et al. (1985, hereafter LLPR). It is important 
to note that the conductivity decreases as the temperature increases. 

For a star with a core temperature ~ 10 8 i^, the important neutrino emission processes 
include electron bremstralung (Flowers and Itoh 1976,1979), neutron-neutron, neutron- 
proton bremstralung and the modified Urea process (e.g. Fridman and Maxwell 1979). The 
rapid cooling processes, e.g. pion condensation (e.g. Maxwell 1977) or direct Urea process 
(Prakash et al 1992) etc, will not be considered here. 

The heat capacities in the crust region come from those of extreme relativistic degen- 
erate electrons and those of nonrelativistic neutrons and ions. The capacities of electrons 
and neutrons are given by Glen and Sutherland (1980). The capacity of ions is given by 
Van Riper (1991). The capacities of the core are mainly due to the relativistic electrons, 



-7- 



the superconducting protons and the superfluid neutrons (Maxwell 1979). The transition 
temperature of normal-superfluid neutrons is given by Takatsuka and Tamagaki (1971). 

We use an empirical formula (Gudmundsson 1983) to relate the surface temperature T s 
and the boundary temperature T b at p& 

T b8 = 1.3(T s 4 6 Mi4)°- 455 , (11) 

where g s is the surface gravity. Eichler and Cheng (1989) have shown that the thermal 
response time between the boundary and the surface is of order of seconds, which is even 
less than our numerical time step (greater than 10 s). The above relation is clearly valid in 
our computation. 

In this paper, we study the thermal response of a glitch which releases heat in a small 
volume inside the star. Although the stellar structure may be spherical symmetric, the heat 
transport is not spherical symmetric as the heat is deposited off center. The pulsar glitch 
may be caused by the sudden transfer of angualr momentum from the more rapidly rotating 
crustal superfluid to the solid crust region or the sudden fracture of the crustal lattice due to 
gravitational stress or magnetic stress. The two heating mechanisms have been studied by 
Alpar et al. (1984a,b) and Ruderman (1976,1991). The energy AE is most likely deposited 
within the density range 10 12 to 2.4 x 10 14 g/cm 3 . AE is of order 10 41 ~ 10 43 ergs (Van Riper, 
Epstein, and Miller 1991). In this paper, we will focus on the case of the Vela pulsar, with 
AE = 10 42 ergs released at ~ 10 l2 g/cm 3 . 

At time t = 0, the heat is deposited in a small volume at r = r g and 9 = 0. r g is 
the radius where the density is p g = I0 12 g/cm 3 (cf. Fig.l). Without lose of generality, we 
choose the hot spot to be centered at 6 — 0, which gets rid of the cp dependence in the 
transport problem. The transport equations in the crust region are solved by explicit finite 
differencing, while the core is taken to be isothermal with a time dependent temperature 
T core determined by the total heat inflow and the core heat capcity. In a typical run, the 
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spatial resolution is taken to be N r x Nq = 100 x 50, with a variable time step determined 
by stability requirements. Numerical convergence of the results have been carefully checked. 
The temperature is cell-centered, while the flux is centered on the cell surface. The inner 
and outer boundary conditions for the flux are given explicitly by 

F r (l,j) = -Ke^~ K - 1 ,3> 1<J<N (12) 

Ar 

at the inner boundary, and 

2 

F r (N r + l,j) = aT*(N r ,j)^ l<J<N e (13) 



" s 

r ob 



at the outer boundary. r b and r s are respectively the radii of the outer boundary of the 
crust and the star, a is Stefan-Boltzmann constant, T s (N r ,j) is the surface temperature of 
the jth angular cell at r b. The initial temperature distribution within the crust is that of 
the equilibrium state of the same star with the initial core temperature. The temperature 



evolution is given by standard finite differencing of Eq. (p~0|) . With a second order scheme 



no extra boundary condition is needed for Eq. ([T0|) as the temperature is cell centered. 



4. Results 

We take the energy to be released at 3 x 10 u g/cm 3 < p < 3 x 10 12 g/cm 3 with a solid 
angle of 2° x 2° centered at 6 = 0°. The initial temperature of the core is taken to be 
T c = 10 8 -fT and the released energy AE = 10 A2 ergs, which is about the energy released 
by the Vela pulsar after a giant glitch. Figure 1 shows the surface temperature versus the 
polar angle at different times. The hot spot is gradually spreading away from 6 = 0° and 
takes about 275 days for the surface temperature to reach the maximum for a UT star. 
The temperature at maximum is about 5 times higher than the background temperature. 
However, the hottest region is concentrated in a solid angle of 2° x 2° centered at 9 = 0°. 
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Figure 2 shows the evolution of the surface luminosity as function of time for three 
different EOSs, i.e. PPS, UT and BPS. The peak luminosity is higher, the time needed 
for reaching the peak is shorter and the relaxation time is also shorter for softer EOS in 
comparison with the stiff one. This is because the soft EOS has a much thinner crust. We 
can see that there are substantial differences among these EOSs. 

Figure 3 compares the luminosity evolutions between the spherical symmetric case and 
non-spherical symmetric case. There are three major differences between these two cases. 
(1) The surface luminosity of the spherical symmetric case reaches the peak almost 5 times 
faster than that of the non-spherical symmetric case. This results from the fact that the 
non-spherical symmetric case has much higher temperature which decreases the conductivity 
a lot. (2) The peak of the total luminosity in the spherical case is lower than that of the 
non-spherical symmetric case. It is because in the non-spherical case more energy comes out 
through the surface instead of heating the core. (3) Most importantly although the changes 
of the total luminosity for spherical case and non-spherical case are very small, the surface 
temperature of the non-spherical symmetric case at a 2° x 2° cap area changes drastically (c.f. 
Fig.l and Fig. 3c of Chong and Cheng 1994 , hereafter CC94) while the surface temperature 
changes very little in the spherical symmetric case. In other words, a hot spot should show 
up on the surface of neutron star 200 to 300 days after the glitch and lasts for about a few 
hundred days if the EOS is soft or intermediate stiff. 

In calculating the evolution of total luminosity for different core temperatures, we find 
that the surface luminosity of the cooler model reaches its maximum earlier. This results 
from the fact that the conductivity is higher for lower temperature. The relative increase 
in the luminosity of the cooler star is higher than that of the hotter star. The net increase 
of the luminosity of the cooler star is a little higher than that of the hotter star. It is 
because the energy of the glitch spreads to stellar surface of the cooler star faster than that 
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of the hotter star. However, the total energy emitted during the thermal afterglow period 
is about the same in these two cases. We have also calculated the thermal evolution inside 
the neutron star at 9 — 0°(cf. Fig.l of Cheng and Li 1997). In comparing with the spherical 
summetric cases (e.g. Figure 2 of CC94), we find that the temperature of the heat pulse is 
much higher in this case because the energy is released in a much smaller volume. The heat 
pulse is propagating outward as well as inward, with a speed slower than that of Figure 2 
of CC94 (where the core temperature is chosen to be 10 6 K, and hence the conductivity is 
much larger). 



5. Discussion 

Based on a set of general relativistic thermal transport and energy balance equations, 
we studied the thermal evolution of a neutron star after a glitch. We find that if the energy 
is released in a compact region, a hot spot can appear on the stellar surface. For a UT 
star with an interior temperature ~ lO 8 -^, although the surface luminosity only increases 
by ~ 10%, the radiation is emitted from a small area with a temperature higher than the 
background temperature by a factor of ~ 5. This results in a periodic hard X-ray pulse 
emission which should stand out clearly from the soft X-ray background. A soft EOS greatly 
enhances this effect and a stiff EOS reduces it: the thermal response to a glitch can provide 
important constraints on the EOS. 

The time for reaching the peak luminosity is long, typically ~ 1 year for a UT star, this 
may make it difficult to relate the hard X-ray pulse to the glitch generating it. However, 
since the energy released by each glitch should be at a different place of the star, comparing 
the relative pulse phase difference and the intensity variation of the hard X-rays observed in 
different epochs can provide evidences for this phenomenon. Together with a detailed spectral 
analysis, the EOS could be deduced. A detailed report on this subject and a comparison 
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with observed data will be presented elsewhere. 



- 12 - 



Figure Captions 

Figure 1. — The thermal profiles of the hot spot at the surface of neutron star of UT 
EOS labelled by the time in days after a heat input of AE = 10 42 ergs induced by a glitch 
released at 3 x 10 n g/cm 3 < p < 3 x 10 12 g/cm 3 in non-spherical symmetric case as a function 
of polar angle. The core temperature is T c = 10 8 iC 

Figure 2. — The evolution curves of the total luminosity for PPS, UT and BPS stars 
denoted by solid, dashed and dot-dashed line respectively. 

Figure 3. — The evolution curves of the total luminosity of UT star in spherical sym- 
metric and non-spherical symmetric cases denoted by dashed and solid line respectively. 
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